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ABSTRACT 

We present a simultaneous periodic and aperiodic timing study of the accreting millisecond X-ray 
pulsar SAX J1808. 4-3658. We analyze five outbursts of the source and for the first time provide a 
full and systematic investigation of the enigmatic phenomenon of the 1 Hz flares observed during the 
final stages of some of the outbursts. We show that links between pulsations and 1 Hz flares might 
exist, and suggest they are related with hydrodynamic disk instabilities that are triggered close to the 
disk-magnetosphere boundary layer when the system is entering the propeller regime. 
Subject headings: stars: individual (SAX J1808.4— 3658 ) — stars: neutron — X-rays: stars 



1. INTRODUCTION 

The low mass X-ray binary transient SAX J1808.4- 

] 3658 (hereafter J1808) was discovered with Beppo-SAX 
in 1996 (in 't Zand et al. 1998) and was found to be an 
accreting millisecond X-ray pulsar (AMXP, Wijnands & 
van der Klis 1998) in 1998 when observed with the Rossi 
X-ray Timing Explorer (RXTE). Since then, four other 
outbursts have been observed with RXTE. The X-ray 

[ lightcurves of the 5 outbursts under consideration are 
remarkably similar in shape and duration. The typical 
outburst duration is several weeks with a recurrence time 

I of ~ 2.5 yr; after 1998, outbursts occurred again in 2000, 
2002, 2005 and 2008. The accretion rate increases steeply 
in the first 2-5 days of the outburst {fast rise), then it 
stays relatively high for a few days (peak), reaching at 
most a few percent of the Eddington rate. After this, 
the X-ray flux has a slow decay lasting ^ 10 days, before 

' entering a fast decay stage in which the flux drops in 3-5 
days. The source then enters a low flux state character- 
ized by 3-5 day flares separated by intervals of very low 
luminosity, the re- flaring state that can last for months, 
followed by quiescence. 
J1808 has shown 401 Hz pulsations during all the 

[ outbursts, at all the luminosities observable by RXTE 
( S 10'^'* ergs~^), even during the re-flares (Hartman et 
al. 2008, 2009a). The re- flaring state was observed to 

; last ~ 60 days (MJD 53550 - 53610) in the 2005 out- 
burst, followed by a low luminosity state approximately 
10 times brighter than quiescence that lasted for another 
^ 60 days (Campana et al. 2008). In the 1998 outburst 
the RXTE observations stopped immediately after the 
beginning of the re-flares, while in the 2000 outburst, 
only the re-flaring state was observed (for ^ 100 days, 
see Wijnands et al. 2001). Thanks to the good sensitivity 
of the AMM-Newton and Swift-XRT satellite, Wijnands 
(2003) and Campana et al. (2008) measured a minimum 
luminosity of ~ 5 x 10'^^ ergs~^ between the flares in the 
re-flaring state in the 2000 and 2005 outburst (assuming 
a distance of 3.5 kpc). In 2002 (Wijnands 2004) and 2008 
(Hartman et al. 2009a) the re-flaring state was observed 
with RXTE for approximately 1 month. 

Campana et al. (2008) interpreted the observed low 
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luminosities as a signature of the onset of the propeller 
regime. The propeller regime is characterized by a Ke- 
plerian velocity in the innermost region of the accretion 
disk that is slower than the rotational velocity of the 
neutron star magnetosphere. In its original formulation 
(lUarionov & Sunyaev 1975) it was proposed to suppress 
the accretion flow onto the neutron star surface. The 
gas, carrying part of the neutron star angular momen- 
tum, was thought to be expelled and spin down the neu- 
tron star. Ghosh & Lamb (1979a, b) proposed that spin 
down and accretion could occur simultaneously and re- 
cent MHD simulations (Romanova et al. 2005, Ustyugova 
et al. 2006) show that two different propeller regimes are 
possible: a strong propeller, characterized by a strong 
outflow of gas, and a weak propeller, with no outflows. 
In both cases a magnetically channeled accretion flow 
onto the neutron star surface is still expected consistent 
with the 401 Hz pulsations observed. 

Wijnands (2004) reported a modulation at a repeti- 
tion frequency of ~ IHz that completely dominates the 
hghtcurve of J1808 in the 2000 and 2002 outbursts. This 
^ 1 Hz modulation appears as sudden intensification of 
the X-ray flux that are obvious in the power spectra and 
sometimes are directly detected in the lightcurve. A re- 
analysis of these data, along with a complete investiga- 
tion of this phenomenon for the other three outbursts 
(1998, 2005 and 2008), is presented in this paper. 

The mechanism of the re-flares and the 1 Hz mod- 
ulation is still unclear, but it might be related to the 
onset of instabilities expected for sources near the pro- 
peller stage (Ustyugova et al. 2006). It has been sug- 
gested that the fast decay and the re-flares are related 
to cooling and re-heating fronts propagating through the 
disk (Dubus, Hameury & Lasota 2001). The heating 
fronts change the accretion disk structure from a neu- 
tral to an ionized state, increasing the viscosity and the 
mass transfer rate through the inner accretion disk. If 
the inner disk structure is influenced by this process, 
the disk-magnetospheric boundary and/or the accretion 
process can be modified as well, possibly producing hy- 
drodynamic instabilities in the accretion flow (Goodson, 
Winglee & Bohm 1997; Spruit & Taam 1993; Bildsten & 
Cutler 1995). 

Whether the 1 Hz modulation is created by such insta- 
bilities is still an open question. J1808 provides a unique 
opportunity to study this, since it shows X-ray pulsations 
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that can be observed simultaneously with the aperiodic 
variability. One of the reasons why the pulsations can 
play an important role in understanding the mechanism 
behind the 1 Hz modulation is the pulse behavior ob- 
served during the 2002 outburst. A drift of ~ 0.2 cycles 
was observed in the pulse phases of the fundamental fre- 
quency (but not in the first overtone), starting and end- 
ing in coincidence with the beginning and the end of the 
fast decay (Burderi et al. 2006). The pulse phase starts to 
drift in coincidence with the beginning of the fast decay 
and ends when the re-flares appear. The interpretation 
of this drift is controversial. 

Burderi et al. (2006) concluded that the phase drifts 
appear in coincidence with the onset of instabilities in- 
duced by accretion of matter onto a weakly magnetized 
star, such as motions of the hot spot on the neutron star 
surface. Hartman et al. (2008) concluded that the ob- 
served phase drift might have been due to a motion of 
the hot spot toward the magnetic pole as the inner ac- 
cretion disk recedes with decreasing luminosity (and thus 
decreasing mass accretion rate). 

In this paper, we present the first comprehensive anal- 
ysis of the 1 Hz modulation and its relation to the re- 
flares and 401 Hz pulsations in all outbursts of J1808. 
We discuss possible explanations for the onset of the 1 
Hz modulation and possible reasons why it has been ob- 
served only in J 1808 until now. We suggest the onset 
of accretion flow instabilities when J1808 enters the pro- 
peller stage as the origin of the 1 Hz modulation. 

2. X-RAY OBSERVATIONS AND DATA REDUCTION 
2.1. RXTE observations 

We reduced all the pointed observations with the 
RXTE satellite's Proportional Counter Array (PCA, Ja- 
hoda et al. 2006) that cover the outbursts of J1808. 

The aperiodic timing analysis was done using 
GoodXenon data with a time resolution of 2^^"s and 
Event data with a time resolution of 2~^^s. The data 
were binned into 1/8192 s bins including all 256 energy 
channels. We performed fast Fourier transforms (FFTs) 
of 128 s data segments, fixing the frequency resolution 
and the lowest available frequency to 1/128 Hz; the high- 
est available Fourier frequency (Nyquist frequency) was 
4096 Hz. No background subtraction or dead-time cor- 
rection was made prior to the FFTs. The Poisson level 
was subtracted from the resulting power spectra. Follow- 
ing Klein- Wolt et al. (2004) we first estimated the Pois- 
son noise using the Zhang et al. (1995) formula and then 
(after inspecting the high frequency range and finding 
no unexpected features) shifted it to match the level be- 
tween 3000-4000 Hz, where no intrinsic power should 
be present, but only counting statistics noise. Then we 
normalized the power spectra using the rms normaliza- 
tion (van der Klis 1995). In this normalization, the in- 
tegral over the power spectrum is equal to the fractional 
rms amplitude squared. The power density units are 
(rms/mean)^Hz~^ and the fractional rms amplitude in 
one specific band is: 

1/2 



TABLE 1 

RXTE AND 5m/t-XTE observations analyzed for each 
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RXTE 










1998 


50914.8 


50939.6 


160 


30411 


2000 


51564.1 


51604.6 


130 


40035 


2002 


52562.1 


52604.8 


664 


70080 


2005 


53522.7 


53587.9 


267 


91056,91048 


2008 


54731.9 


54776.1 


236 


93027 


SWIFT 










2008 


54756.6 


54778.3 


15 


0003003435-00030034044 
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The errors on the fractional rms are calculated by us- 
ing the dispersion of points in the data. We consider a 



measurement as a non-detection in a specific band, when 
the ratio between the fractional rms and its standard de- 
viation is smaller than 3. In this case we quote upper 
limits at the 98% confidence level. 

The periodic 401 Hz pulsations were analyzed by con- 
structing pulse profiles folding chunks of lightcurvc in 
profiles of TV = 32 bins, with the cphcmcris of J1808 
provided by Hartman et al. (2009a). In this folding 
process we used the TEMPO pulsar timing program (v 
11.005) to generate a scries of polynomial expansions of 
the ephemeris that predict the baryccntered phase of 
each photon detected. The length of each data chunk 
was chosen according to the length of each RXTE obser- 
vation (Obs-Id). 

We then split each pulse profile into a fundamental (at 
v) and a first overtone (at 2v) using standard fits. For 
a detailed description of the method used for measuring 
the pulse time of arrivals (TOAs) we refer to Patruno 
et al. (2009). A set of pulse phase residuals was then 
obtained by subtracting a Keplerian circular orbit and 
constant pulse frequency model, using the ephemeris of 
Hartman et al. (2009a). 

2.2. Swift-XRT observations 

During the 2005 and 2008 outbursts, the Swift X-ray 
Telescope {Swiftr-XKT) observed the re-flaring state of 
J1808. We refer to Campana et al. (2008) for the study 
of the Swift>-XKT 2005 observations. Here we consider 
10 pointed observations (see Table 1) covering a total 
of ^ 15 ks that were taken during the 2008 outburst 
(MJD 54757-54778) and were reduced by using the XRT 
pipeline (v. 0.12.0). Each observation lasted between 1 
and 3 ks. The data were collected in photon counting 
(PC) mode, except for observation 00030034041 which 
was taken in windowed timing (WT) mode. 

We extracted source and background events for each 
observation using circular regions with radii of 20 arc- 
seconds, and extracting photons with energies between 
2-10 keV and 0.5-10 keV. 

3. RESULTS 

3.1. The X-ray lightcurves and the re-flaring state 

In the 2008 outburst rc-flaring state, the lightcurvc 
reached very faint luminosities, with a large portion of 
the re-flares below the sensitivity limit of RXTE, but 
above the detection threshold of Swift-XRT. The flux 
reached a minimum level of 2.0 x 10^^'' erg s^^ cm~^ in 
the 0.5-10 keV band, which corresponds to a luminosity 



Re-Flaring state, 1 Hz QPO and pulsations 



3 



of ^ 3 X 10'^^ erg s ^ at 3.5 kpc. This luminosity is of the 
same order of magnitude as that observed by Wijnands 
(2003) and Campana et al. (2008) during the 2000 and 
2005 re-flares. 

A property of all the re-flares is the periodicity on time 
scales of a few days, creating the "bumps" observed in the 
lightcurve (see Fig. 1). The fast decay and the re-flares in 
the 2002 and 2005 show a similar duration. The 2008 out- 
burst also shows comparable timescales, although with 
much higher uncertainty. As noted in § 1, in '98 the ob- 
servations stopped too early and the '00 ones started too 
late to allow a similar comparison. 

The precise determination of the re-flare periods suf- 
fers of biases due to occasional poor sampling and to 
the sensitivity limit of RXTE. Therefore, although the 
data certainly allow this, there is no significant evidence 
that differences in fast decay time scales between the 3 
outbursts have an effect on the subsequent re-flare peri- 
odicity time scales. 

In the 2008 outburst it is hard to calculate the re- 
flare time scale because the flux is below the sensitivity 
limit of RXTE in the majority of the observations. The 
Swift sampling was approximately 1 observation every 
2 days, too long compared to the fast decay time scale 
of 3 days to unambiguously exclude shorter time scales. 
However, the fast decay and the re-flare time scales are 
again compatible with being the same. 

In 2002 and 2005 the re-flares' peak luminosity tends 
to decrease with time (see Fig. 1). This is not observed 
in the 2000 (see Wijnands et al. 2001) and 2008 re-flares, 
which show an erratic change of peak luminosities. In 
the re-flaring state, the luminosity can change by ~ 3 
orders of magnitude on time scales of ^1.5 days (see 
Fig. 1 at MJD-^ 54756 and Wijnands 2003 for a similar 
observation in the 2000 outburst). 

3.2. The fast decays and the pulse phase drifts 

Close to the end of the 2002 outburst, the pulse phase 
of the fundamental was observed to drift by 0.2 cycles in 
just a few days, in coincidence with the beginning and 
the end of the fast decay (Burderi et al. 2006). Hartman 
et al. (2008) showed that a similar drift was present in 
2005. Here we show that both phase drifts start and end 
in coincidence with (or very close to) the beginning and 
the end of the fast decay. 

In Fig. 2 we plot the pulse timing residuals of the fun- 
damental and the 2-10 keV X-ray flux in mCrab. In the 
2002 outburst the pulse phase drift starts exactly when 
the fast decay begins. In the 2005 outburst we see the 
phase drift beginning close to the beginning of the decay, 
although given the large scatter in the pulse phases it is 
difficult to define the exact moment when the phases be- 
gin to drift. Both the 2002 and 2005 pulse phase drifts 
end in coincidence with the end of the fast decay. Both 
the 2002 and 2005 pulse phases drift by ~ 0.2 cycles. 

The slope of the fast decay is consistent with being 
the same in both outbursts (see also Fig. 3). The pulse 
phases in both outbursts during this time drift at the 
same speed of ^ 0.07 cycle/day. The pulse phase behav- 
ior in the 2002 and 2005 fast decay is therefore consistent 
with being identical. 

In the 2000 outburst a similar test cannot be performed 
since RXTE vsAssed the decaying portion of the outburst, 
and the observations covered just the re-flaring state. In 
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FlK. 1. - Outburst lightcurvcs in the 2-10 kcV energy band. 
The count rate was iiomialized to the Crab intensity measurements 
nearest in time (Kuulkcrs et al. 1994) and in the same PCA gain 
epoch (e.g., van Straaten et al. 2003). Each black circle and triangle 
corresponds to one RXTE andSu)i/t-XRT observation, respectively 
(see Table 1 for a complete list of the observations used) . The open 
squares indicate observations when the 1 Hz QPO was detected, 
while the open circles indicate the RXTE observations with an X- 
ray flux below the detection threshold of ~ 1 mCrab. The vertical 
dashed lines divide the lightcurve into intervals beginning with the 
fast decay, and separating the re-flares. The duration of eax;h re- 
flarc is indicated in days and is similar to that of the fast decay 
(also indicated). In the 2008 outburst only few observations were 
above the detection threshold and all the vertical lines refer to the 
fast decay time scale. 

the 1998 and 2008 outbursts no clear phase drift is ob- 
served. We refer to Hartman et al. (2008, 2009a) for 
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Fic;. 2. — Pulse phase residuals (upper panels) calculated for the fundamental frequency alone. The bottom panels show the X-ray 
lightcurvc (in logarithmic scale). The dashed lines mark the beginning of the fast decay in the lightcurves. At this point the pulse phases 
start to drift for ~ 0.2 spin cycles. The solid black lines mark the time of the first appearance of the 1 Hz QPO. Note that in both outbursts 
it appears at the end of the phase drift close to (but before) the minimum flux level in the X-ray lightcurve. 




Fig. 3. — X-ray lightcurve of four outbursts, plotted up to the end of the fast decay stage. The curves are aligned to the beginning of 
the fast decay. The 2008 outburst has the shortest decay time. The 2008 slow decay, peak and fast rise stages are very similar to those of 
the 2005 outburst. The end of the 1998 outburst fast decay is limited by a non detection (open circle); it is the outburst that shows the 
dimmest luminosities before the beginning of the re-flares. 



a detailed discussion of the coherent timing analysis of 
those two outbursts. 

3.3. QPO parameters and flux 

We will not in this paper provide a complete descrip- 
tion of all the aperiodic variability observed in the 5 out- 
bursts. For a description of the aperiodic timing fea- 
tures observed in AMXPs, we refer to van Straatcn ct al. 
(2005) and to the review of Wijnands (2006). Here we fo- 
cus our attention on the 1 Hz modulation and give a brief 
description of the power spectra observed from the fast 
decay stage on, in order to provide context. A descrip- 
tion of our quantitative analysis of the IHz modulation 



follows in the next sections. 

We examined the entire data set of all 5 outbursts for 
evidence of the 1 Hz modulation. Close to the end of the 
2002 and 2005 pulse phase drifts, when the fast decay 
state is almost over, the X-ray lightcurve clearly shows 
this strong modulation with a repetition frequency of ~ 1 
Hz (Fig. 4) . It shows up as a strong quasi periodic oscil- 
lation (QPO) peak around 1 Hz in the power spectrum 
(Fig. 5, see also Wijnands 2004 for a similar plot for the 
2002 outburst). Later during the re-flares of these out- 
bursts the modulation occasionally recurs, as discussed 
in more detail in § 3.4. The 1 Hz modulation also ap- 
pears and disappears sporadically during the re-flares of 
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the 2000 but was not detected in the 1998 and 2008 out- 
bursts (cf. § 1). 

We calculated the fractional rms amplitude^ of the 1 
Hz QPO by integrating the power in the 9.95 Hz wide 
band 0.05 — 10 Hz. This frequency band was chosen be- 
cause it is here that the power is observed in the large 
majority of the observations. Some power above 10 Hz 
is observed in some cases as an extended tail of the 1 Hz 
QPO (Fig. 5), but selecting an upper frequency of 20 Hz 
does not change the results significantly. Clearly, using 
this method we sum up all the power from the different 
components in the power spectrum: power of the 1 Hz 
QPO, its harmonics, and the broad band noise. So every- 
thing in the range 0.05-10 Hz is included, and therefore 
this method is not optimal to measure the power gener- 
ated by one single component. However, the rms calcu- 
lated in this way is empirical and independent from any 
model used to fit the data. The fractional rms amplitude 
in the 0.05-10 Hz band was in the range 10—125% in all 
observations where the modulation was detected. 

The QPO is usually quite broad, and its shape cannot 
be satisfactorily fitted by single or multiple Lorentzians, 
as it has a sharp fall-off in power at lower frequencies. 
In some observations a second harmonic peak is visible 
(Fig. 5). We model these features with Gaussians of the 
form: 

P{vi)=A - e^^[{vi-Bf/C^] (2) 

where Vi are the Fourier frequencies, A is a normalization 
factor, B is the ccntroid frequency (i^o or 2z/o) and C is 

related to the FWHM through: FWHM = 2C • [ln(2)]^/^ 
The typical value for the FWHM of the two Gaussians is 
^ IHz for all the three outbursts. The Gaussian that we 
use to fit the 1 Hz modulation usually has a quality factor 
Q = i/q/FWHM < 2. When Q < 1 the power spectrum 
does not show a clear peak, but sometimes it does show 
a break in the noise (see for example the bottom panel 
in Fig. 5). Here we call the 1 Hz modulation "QPO" 
regardless of its quality factor. When Q < 1 wc did 
not try to fit the power spectra with gaussians because 
the centroid frequency is ill defined. When the 1 Hz 
QPO is fitted by 2 Gaussians vq is defined as being the 
frequency of the Gaussian with the highest peak power 
in the power spectrum. With this choice is consistent 
with being always the fundamental frequency of the 1 
Hz modulation (varying between 0.8 and 1.6 Hz), as the 
harmonic peak at 2v(j is always lower in maximum power. 

We calculated the fractional rms amplitude and cen- 
troid frequency of the 1 Hz QPO and the X-ray flux 
for the 2000, 2002 and 2005 outbursts. The fractional 
rms amplitude as observed in 2002 and 2005 is shown 
in Fig. 6 (middle panel). The fractional rms amplitude 
shows abrupt changes over time and does not follow a 
clear correlation with flux. In the 2000 re-flaring state, 
the 1 Hz modulation was detected in 13 observations out 
of 46, close to the peaks of the re-flares. 

We plotted together all the points for the three out- 
bursts in Fig. 7. There is a clear increase of the fre- 

^ Fractional rms amplitude is standard deviation divided by 
mean flux. Values of rms larger than 100% indicate that the 1 
Hz flares have a large amplitude and a short duty cycle For a 
light curve composed exclusively of square flares with duty cycle 

/, rms = J , arbitrarily large for / — ► 



quency with flux and and an anti-correlation between 
frequency and QPO fractional rms amplitude. Both re- 
lations have considerable scatter. We performed a rank 
correlation test on the 1 Hz QPO frequency vs. frac- 
tional rms anticorrelation and we found a Spearman co- 
efficient oi p = —0.38 with a probability of 0.1% of the 
null hypothesis (no correlation in the data) being true. 
A similar test for the frequency vs. flux correlation gives 
p = 0.75 with a probability of less than 0.01% for the 
null hypothesis. 

We explored the fractional rms dependence of the 1 
Hz QPO on the X-ray flux using a similar plot, and we 
found no clear dependence. 

We then calculated the upper limits (quoted at the 
98% confidence level) for all the observations after the 
end of the fast decay where the QPO was not detected. 
The upper limits were calculated per observation. In 
the 1998 data two complications occur: the observations 
ended immediately after the beginning of the re-flares, 
and the fast decay reached very faint fluxes (< 1 mCrab) 
at its minimum. Only in one observation J1808 was de- 
tected; then the upper limit was 55% rms. A 1 Hz mod- 
ulation therefore cannot be completely excluded for this 
outburst. In 2000, 2002, 2005 and 2008 the most con- 
straining upper limits are 15%, 9%, 19% and 19% rms, 
respectively. 

Across the 2000, 2002 and 2005 outbursts, the 1 Hz 
QPO was observed in a rather narrow range of luminosi- 
ties during the re-flares (2-15 mCrab, 2-10 keV), in con- 
trast with the larger range of luminosities covered by the 
outbursts (~ 1-80 mCrab). The SwiftrXRT and XMM- 
Ncwton observations had insufficient time resolution or 
an insufficient number of counts to probe the presence of 
the 1 Hz QPO below 1 mCrab. 

3.4. The appearance of the 1 Hz QPO 

In the 2002 outburst, the IHz QPO appears immedi- 
ately after the pulse phase drift is complete, but with the 
flux still decreasing in the fast decay stage. In 2005 the 1 
Hz QPO appears in a similar position (see black vertical 
lines in Fig. 2). In both the 2002 and 2005 outbursts, 
data gaps prevent the observation of the exact moment 
when the 1 Hz QPO appears. Taking account of the 
gaps, the flux level of the first 1 Hz QPO appearance is 
consistent between the two outbursts. 

To track the appearance of the 1 Hz QPO on time 
scales as short as a few hundred seconds, we investigated 
the observation available in which the 1 Hz QPO first 
appears in the 2005 outburst. The observation (Obs-Id 
91056-01-04-03) shows two chunks of data separated by 
a gap of 5000 s (chunk A with a length of 3200 s and 
chunk B with length 3500s) that we analyzed separately. 

In chunk A the power spectrum has already changed 
with respect to the earlier observations: the power at 
higher frequencies has disappeared, and there is only 
power in the range 0.05-10 Hz. Probably we are witness- 
ing the onset of the 1 Hz modulation. The average power 
spectrum shows a fractional amplitude of (21 ± 1)% rms 
in the 0.05 10 Hz band. In chunk B we find a 1 Hz QPO 
clearly present with a fractional amplitude of (68 ±4)% 
rms. 

We then calculated power spectra from 256 s data seg- 
ments and determined the fractional rms amplitude in 
the 0.05-10 Hz band, (see Fig. 8). Significant power 
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Fig. 4. — 1 Hz modulation of the X-ray lightcurve (background subtra<;ted) , as observed in a re-flare of the 2005 outburst. The figure 
shows a 15 s chunk of lightcurve with a time resolution of 0.1 seconds. In this observation the 1 Hz modulation has the highest measured 
fractional rms amplitude (125% rms, Obsld:91418-01-02-05). 



is detected in three segments during the observations of 
chunk A, meaning that some low level activity is already 
present in the same frequency range where the QPO will 
later appear. In chunk B power is nearly always detected; 
the onset of the 1 Hz QPO must have occurred within the 
5000 s gap. At MJD 53540.48 the amplitude shows an 
abrupt decrease from ^ 60% rms down to < 30% (98% 
confidence level) on a time scale of 256 s. 

We also calculated the power in a number of frequency 
bands in the range 1/128 to 256 Hz for all power spectra 
of chunk A. The power is always consistent with zero at 
the three sigma level in all the frequency bands except 
for the 0.05-10 Hz band in the three segments mentioned 
above. 

3.5. Energy dependence of the 1 Hz QPO 

In Fig. 9 wc show the energy dependence of the 1 Hz 
QPO. The points in the figure refer to the observation 
(one per outburst) for which the fractional rms of the 1 
Hz QPO in the 2-60 keV energy band was the highest 
(110.5% rms in 2000, Obsid 40035-01-04-01, total expo- 
sure 7 ks; 117% rms in 2002, ObsId 70080-03-15-00, total 
exposure 2 ks; 125% rms in 2005, ObsId 91418-01-02-05, 
total exposure 1.2 ks). In these observations the char- 
acteristic frequencies of the 1 Hz QPO and its overtone 
form clear peaks in the power spectrum and Q > 2. 

The energy spectrum is hard, rising by a factor of 1.5- 
1.7 between 2 and 17 keV. Upper limits (98 % confidence 
level) in the 17-60 keV band were: 114% rms for 2000, 
157% rms for 2002 and 155% rms for 2005. 

3.6. Jumps in 401 Hz pulse phases related to strength of 

1 Hz QPO 

During the first re-flare of the 2002 outburst, the 1 Hz 
QPO becomes very broad and then disappears. When 



the QPO becomes very broad and its fractional rms am- 
plitude suddenly drops from ~ 50% down to 10%, a 
jump of 0.2 cycles is observed in the pulse phases of the 
fundamental (jump A, see Fig. 6). Two similar pulse 
phase jumps of 0.1 cycles are observed (jump B and C) 
when the QPO is not detected, with rms amplitude up- 
per limits of ~ 30 and ~ 10% rms (98% confidence level). 
We call these phase changes "jumps" as opposed to the 
"drifts" observed during the fast decay stage (§ 3.2), since 
the 0.1 — 0.2 cycle jumps are sudden, occuring on time 
scales of a few hours or less, while the 0.2 cycle drifts 
take 4 — 5 days. 

In the 2005 outburst we see something very similar, al- 
though the re-flares have a sampling that is much worse: 
the observations are rather sparse and are separated by 
at least 1 day. In both the 2002 and 2005 outbursts, the 
pulse phase in the jumps is close to the pulse phase prior 
to the beginning of the fast decay stage (see § 3.2). 

In the 2000 outburst the observations arc rather sparse, 
with gaps of several days between observations. However 
also in this case when the QPO rms amplitude is low or 
the QPO is undetectable, the pulse phases are observed 
to jump by ~ 0.1-0.2 cycles, consistent with the 2002 
and 2005 pulse phase behavior. 

This phenomenology relating the 1 Hz QPO rms am- 
plitude and the pulse phases might suggest a link be- 
tween the presence of the 1 Hz QPO and the accretion 
flow onto the neutron star surface. However, due to the 
sparse sampling of the observations and the small num- 
ber of jumps observed, this link remains to be confirmed 
in future observations. 

3.7. Dependence of 401 Hz pulsation on the 1 Hz phase 

The 1 Hz modulation might affect the magnetic chan- 
neling during the re-fiaring state. Possible differences in 
the pulse properties can be unveiled by analyzing the pul- 
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Fig. 5. — Power spectra (power X frequency) of the rc-flares of 
the 2005 outburst. The four plots show different manifestations 
of the 1 Hz modulation in the power spectrum. The upper panel 
(Obsld:91418-01-02-0r3, MJD~ 53550.1) shows a strong coherent 1 
Hz QPO and its overtone clearly separated. They blend together in 
the second panel (Obsid: 91418-01-01-00, MJD~ 53542.5) where 
some power is also observed at frequencies higher than 10 Hz. The 
third panel (ObsId: 91056-01-04^01, MJD~ 53541.5) shows a QPO 
with a steep cutoff at frequencies larger than ~ 10 Hz. The bottom 
panel (ObsId: 91418-01-03-04, MJD~ 53556.7) shows still power, 
but as an incoherent feature without a clear peak. 

sations coming from different phases of the 1 Hz modula- 
tion. Since the modulation has a frequency that changes 
with time, it is hard to unambiguously define its phase. 

We therefore used a crude but simple approach, defin- 
ing "peaks" and "valleys" of the 1 Hz modulation in the 



lightcurvc relative to the average flux. We first calculated 
the average X-ray flux per observation. Then wc split the 
X-ray lightcurve into an upper half and a lower half con- 
taining peaks and valleys, respectively. We checked that 
this method was properly dividing the lightcurve by a 
visual inspection. We then folded all the peaks and the 
valleys into two different 401 Hz pulse profiles at the spin 
frequency of the neutron star (see Hartman ct al. 2008 
for the ephemeris used to fold the data), and repeated 
the procedure for each observation where the 1 Hz QPO 
was detected. 

The pulse phases are always consistent with being the 
same for peaks and valleys, within the pulse phase un- 
certainty of down to 0.01 cycles. The absolute pulse am- 
plitude is higher in the peaks than in the valleys but by 
a smaller fraction than the mean flux, so that the pulse 
fractional amplitude is always higher in the valleys by 
a factor varying between 1.2 and 6 (with a typical error 
on the amplitude ratios of ~ 0.1). The pulse amplitudes 
of peaks and valleys follow a similar evolution in time: 
when the pulse amplitude of peaks rises (or drops) in an 
observation, the same is seen for the pulse amplitude of 
valleys. 

The valleys fractional rms amplitudes span a range be- 
tween a few % rms up to more than 10 % rms, with val- 
ues much larger than the average pulsed fractions of the 
persistent pulsation observed during the outbursts before 
the reflares (for a discussion of the pulse fractional ampli- 
tudes during the 1998, 2000, 2002 and 2005 outbursts see 
Hartman et al. 2008 and for the 2008 outburst see Hart- 
man et al. 2009a). The opposite is true for the pulsed 
fractions of peaks: they are on average smaller than the 
observed fractional amplitudes during the outbursts. 

4. DISCUSSION 

Our analysis of the five outbursts of SAX J1808.4-3658 

has shown that the fast decay, the re-flares and the 1 
Hz modulation can be related phenomena in the 2002 
and 2005 outbursts (and possibly in the 2000 one) which 
might also affect the 401 Hz pulsations. The drifts in 
pulse phase and the re-flares observed on time scales of 
days might be related with some change in the disk struc- 
ture which is connected to the fast decay, which takes 
place on a similar time scale. The fast decay presumably 
is either a thermal or a viscous one, and additionally 
depends on the lengthscale on which significant changes 
occur within the disk, longer length scales corresponding 
to longer time scales. The 1 Hz QPO properties had pre- 
viously only been very briefly discussed in the literature, 
and left the mechanism responsible to be identified. We 
therefore studied the 1 Hz modulation properties in de- 
tail and related them to the behavior of the X-ray flux 
and the X-ray pulsations. Any suitable model for the 1 
Hz modulation has to explain the following key proper- 
ties that we reported in § 3: 

• In 2002 and 2005, it appears at the end of the fast 
decay, after the pulse phases have drifted by 0.2 
cycles (Fig 2) 

• It appears at a flux level that is consistent in the 
2002 and 2005 outbursts (§ 3.2) 

• It appears sporadically only in a narrow range of 
luminosities (^15mCrab, § 3.3). 
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• It is very coherent in some observations while in 
some others it is a broad incoherent feature (Fig. 4 

and 5) 

• In a few cases, its amphtude might be connected 
with the 401 Hz pulse phases (§ 3.6, Fig. 6). 

• Its fractional rms amplitude can be as high as 
125% and then 10% for the same X-ray luminos- 
ity (Fig. 6) 

• Its amplitude is energy dependent, rising with en- 
ergy up to ~ 17keV (Fig. 9) 

• Its rms amplitude is very high, up to 170% in the 
10-16 keV band (Fig. 9) 

• Its centroid frequency is quite stable and slightly 
increases with flux (Fig. 7) 

• Its centroid frequency decreases with increasing 
fractional rms amplitude (Fig. 7) 



Flares that might be similar to the J 1808 re-flares 

arc seen in other low mass X-ray binaries (LMXBs): 
the black hole candidate XTE J1650-500 (Tomsick et 
al. 2003), and the neutron star transient SAX J1750.8- 
2900 (Linares et al. 2008). This behavior therefore may 
not be unique to J1808, although this is the only AMXP 
in which such re-flares have been observed. 

However, the 1 Hz modulation does appear to be nearly 
unique to J1808. The most similar QPOs, in terms of fre- 
quency and high fractional amplitude, are found in the 
systems that show Type II bursts. The Rapid Burster 
(Lewin et al. 1976; Hoffman, MarshaU & Lewin 1978; 
Marshall et al. 1979) has QPOs in the range 0.04-4 Hz 
in the persistent emission, and from 2-5 Hz in during its 
Type II bursts (Tawara et al. 1982; Lewin 1987; Stella 
et al. 1988; Dotani et al. 1990; Lubin et al. 1991, 1992; 
Rutledge et al. 1995). Similar frequencies (0.04-0.4 Hz) 
have also been reported in the aftermath of Type II 
bursts from the bursting pulsar GRO J1744-28 (Kom- 
mers et al. 1997). This similarity may point to a magne- 
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does in any case seem to rule them out. 

There arc various oscillatory modes of the neutron star 
surface layers that could lead to periodic brightness vari- 
ations. The most likely candidates would be an ocean 
g-mode (a vibration driven by thermal buoyancy Mc- 
Dermott & Taam 1987; McDermott, van Horn & Hansen 
1988), although to obtain a 1 Hz frequency would neces- 
sitate a very high order harmonic(Bildsten, Ushomirsky 
& Cutler 1996). The dependence of amplitude on energy 
could be explained by a mode model, but modes cannot 
explain the extremely high fractional amplitudes of the 1 
Hz QPO (Piro & Bildsten 2006). In addition, if this was 
the right mechanism, then it should also be observed in 
other stars since the triggering conditions should not be 
unique to J 1808. 

The matter accreted onto the surface of the neutron 
star may burn stably, unstably (generating X-ray bmsts), 
or in a marginally stable fashion (for a review sec Bild- 
sten 1998). The time scale for the quasi-periodic vari- 
ations associated with marginally stable burning is set 
by the accretion time scale and the thermal time scale. 
For hydrogen ignition, the time scale for marginally sta- 
ble nuclear burning, computed as in Heger, Gumming 
& Woosley (2007) using the appropriate accretion and 
thermal time scales, is clearly too slow (~ 10 minutes). 
For helium ignition, the time scales are shorter but still 
too slow 100 s), thus ruling out the marginally stable 
nuclear burning scenario. In the following sections we 
therefore focus on the accretion flow as the most likely 
mechanism. 



4.1. Accretion onto a magnetized neutron star 

In order to understand what might be causing the 1 Hz 
QPOs, some general background on magnetically chan- 
neled accretion will be useful. To sustain channeled ac- 
cretion at the maximum accretion rate of a few percent 
Eddington, the magnetic field for J1808 must be > 4x10'' 
G (Hartman et al. 2008). The upper limit on the field, 
determined from timing, assuming that the spin down 
comes from magnetic dipole radiation from a rotation 
powered pulsar, is 1.5 x 10^ G (Hartman et al. 2008, 
2009a). 

In a discussion of magnetized accretion, reference is 
often made to the corotation radius rc, the radius at 
which matter in a Keplerian orbit would have the same 
angular velocity as the star. 



tosphere/disk mechanism (S 4.2). 

Oscillations between 0.58 and 2.44 Hz were also re- 
ported in three dipping sources, although with an energy 
independent rms amplitude of below 12% (Homan et al. 
1999; Jonker, van der Klis & Wijnands 1999; Jonker et 
al. 2000). Disk warping or shadowing was suggested as 
the most likely mechanism for these QPOs (see § 4.2.1). 

Alternatively, the 1 Hz QPO may be caused either by 
the accretion flow (including possible occultation phe- 
nomena), or by some process that occurs after matter 
arrives on the surface. If the latter, it is very hard to 
understand why it is only seen in J1808. A brief con- 
sideration of the two most plausible surface mechanisms, 
global oscillations or marginally stable nuclear burning, 
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where Ug is the spin frequency of the neutron star and 

M is its mass. For J1808, assuming a mass of 1.4M0, 
rc ^ 31 km. The other radius of relevance is the mag- 
notospheric radius r„i: the radius at which the magnetic 
field becomes dynamically important in controlling the 
inflow of matter. For spherically symmetric accretion, 
one can estimate by setting the magnetic pressure 
equal to the ram pressure of free fall (Lamb, Pethick & 
Pines 1973): 



10 



Patruno et al. 



'7.8 



" B ' 


4/7 


■ R ' 


12/7 


■ M ' 


-1/7 


108 G 




10 km 




1.4 Mq_ 





M 



-2/7 



km 



(4) 



where we have assumed a dipolc field B ^ f^/r^-, /i being 
the magnetic moment. R is the stellar radius and M 
the accretion rate. In the case where accretion occurs 
from a disk this expression will be slightly modified by 
the rotational energy of the disk (Spruit & Taam 1993). 
Rotation of the central star can also afTcct the location 
of rm (Lovelace, Romanova & Bisnovatyi-Kogan 1999). 

This simple order of magnitude estimate yields 18 km 
for SAX J 1808 at the peak of the outburst (assuming 

M = 5%MEdd and B = 10® G). This increases as the 
accretion rate falls, becoming comparable to Vc once the 
accretion rate drops to 1% of the Eddington rate. The 
precise value of will depend on details of inner disk 
physics as discussed in Psaltis & Chakrabarty (1999). 

When r-m < Vc accretion should proceed without dif- 
ficulty, with the magnetic field channeling material out 
of the disk and onto the magnetic poles (Pringle & Rees 
1972). Once > rc, however, the situation becomes 
more complex, and the system is said to be in the so- 
called "propeller" regime. Initially it was thought that 
for r-m > Tc accretion would cease, with matter being 
expelled from the system (Illarionov & Sunyaev 1975). 
Further study by Spruit & Taam (1993) showed that 
actually has to exceed Vc by a reasonable margin for ma- 
terial to be expelled. Steady accretion is in fact possible 
when Vm > rc, even though the neutron star should spin 
down. In this stage the inner edge of the disk stabilizes 
near Vc and the density at the inner disk rises allowing 
angular momentum to be transferred outwards (Spruit 
& Taam 1993; Rappaport, Fregeau & Spruit 2004). This 
type of disk structure predicts spin-down without requir- 
ing penetration of the disk by the magnetic field far be- 
yond rc, as was proposed in early works on the topic 
(Ghosh & Lamb 1979b). In this paper we use the term 
"propeller" to define the condition in which r.^ > re- 
in the next Section we will examine in turn the various 
mechanisms that might be responsible for the observed 
1 Hz QPO. 

4.2. Candidate Mechanisms: disk/magnetosphere 
instabilities 

There are numerous ways of obtaining variability from 
the disk and its interactions with the stellar magneto- 
sphere. Many have been explored, however, as a means 
of explaining the kHz QPOs - and so have frequencies 
that would be very far off 1 Hz. For this reason we will 
neglect many of the mechanisms for disk variability that 
have been discussed in the literature (see van der Klis 
2006 for an extensive review of these mechanisms) and 
focus on those that might have frequencies in the right 
range. 

4.2.1. Disk obscuration 

The dipping QPOs have fractional amplitudes of ~ 
10% (Homan et al. 1999; Jonker, van der Khs & Wij- 
nands 1999; Jonker et al. 2000). In the dipping sources 



the energy dependence of the QPO amplitude was flat, 
supporting the shadowing hypothesis. In J1808, there 
is a clear energy dependence and the 1 Hz QPO am- 
plitude is much higher than that seen in the dippers. 
Furthermore, there is no evidence for dipping in the X- 
ray lightcurve of J1808. Several recent studies ( Cackett 
et al. 2009, Papitto et al. 2009, Ibragimov & Poutanen 
2009, Deloye et al. 2008) suggest an inclination for J1808 
of around 60° which is too small to produce dipping. The 
mechanism proposed for this QPO was a partial obscu- 
ration of the neutron star surface via a blob of plasma 
in the disk, orbiting at a keplerian frequency of ^ 1 Hz. 
The disk radius corresponding to a keplerian frequency 
of 1 Hz is r ^ 1700 km. Considering the observer at 
an inclination of 60°, the blob of plasma needs to be as 
thick as 1000 km to allow partial obscuration of the neu- 
tron star surface. The expected accretion disk thickness 
at that radius is expected to be however orders of mag- 
nitude thinner than this value. By using for example 
Eq. (9) in Rappaport, Fregeau & Spruit (2004), the disk 
has an expected scale height of ~ 70 km. So, the 1 Hz 
QPO in J1808 is unlikely to be explained with a dipping 
mechanism. 

Another possibility is the occurrence of occulations of 
the neutron star surface when the inner edge of the ac- 
cretion disk enters in the line of sight of the observer. 
The 1 Hz QPO appears only at very low flux levels (in 
the range ^ 0.001 — 0.01% Eddington) which has two 
important consequences. First, by using Eq. 4, the in- 
ner edge of the disk is at approximately 56 km when the 
X-ray flux is at its minimum (0.001% Eddington). This 
means that, by assuming an observer inclination of 60°, 
the neutron star surface is obscured for a disk thickness 
of ~ 60 km. This thickness is however much larger than 
the expected disk scale height which is only 0.04 km. 
Moreover, even considering the physical size of the inner 
disk as large as 60 km, the optical depth would drasti- 
cally reduce to values below 1 after a few km in height 
from the disk middle plane. Although the inner part of 
the accretion disk could be optically thin, allowing pho- 
toelectric absorption and thus explaining the hard energy 
dependence of the 1 Hz QPO rms amplitude, it cannot 
explain the very high rms amplitude of the QPO, since 
the optically thin plasma cannot completely obscure the 
neutron star surface. 

4.2.2. Interchange instabilities 

The role of interchange instabilities in admitting mat- 
ter to the magnetosphere is discussed in detail by Arons 
k Lea (1976) and Eisner & Lamb (1977). The mag- 
netic pressure prevents incoming matter from crossing 
magnetic field lines, but if it is energetically favorable 
(due for example to gravity) for material to be "inside" 
the field lines rather than "outside", then interchange 
instabilities can act to move the material inside. This 
is often referred to as the Rayleigh- Taylor instability, a 
term that refers to the instability that occurs when a 
more dense fluid overlies a less dense fluid in a gravita- 
tional fleld. Where a plasma is supported against grav- 
ity by a magnetic field, this is more correctly referred 
to as the Kruskal-Schwarzschild instability (Kruskal & 
Schwarzschild 1954). 

The conditions necessary for the onset of Kruskal- 
Schwarzschild instabilities in the situation where ^ rc 
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have since been studied using MHD simulations (Ro- 
manova, Kulkarni & Lovelace 2008; Kulkarni & Ro- 
manova 2008). For misalignment angles < 30° accre- 
tion can proceed either stably via funnel flows, or un- 
stably via interchange instabilities. In the unstable sit- 
uation matter accretes via a number of "tongues" that 
penetrate the magnctopausc in the equatorial plane. If a 
certain number of tongues dominate, quasi-periodic os- 
cillations can emerge in the light curves. Funnel flows 
can co-exist with accretion by tongues (Kulkarni & Ro- 
manova 2008), although their presence should reduce the 
amplitude of the persistent pulsations by reducing the 
azimuthal asymmetry (Romanova, Kulkarni & Lovelace 
2008). This might be consistent with the continued pres- 
ence of accretion-powered pulsations, with an amplitude 
that depends on the phase of the 1 Hz QPO. 

In the cases stiidied by Romanova, Kulkarni & 
Lovelace (2008) interchange instabilities set in above a 
critical accretion rate, making them unlikely as a cause 
for the 1 Hz QPO (which appears only below a criti- 
cal rate). Interchange instabilities cannot be completely 
ruled out, however. In the standard case studied by 
Romanova, Kulkarni & Lovelace (2008) magnetic pres- 
sure and gravity dominate the force equations. When 
fm ^ fc, however, magnetic pressure equals gravity. At 
this point other terms start to dominate the force equa- 
tions and the character of the interchange instability will 
change (Baan 1977; Spruit & Taam 1993). Baan (1977) 
showed that sporadic penetration of the magnetosphere 
is possible in this regime. However detailed numerical 
simulations of the type performed at higher accretion 
rates have not been done, and the effect on funnel flows 
(and hence the amplitude of the accretion-powered pul- 
sations) is not known. Without further study, periodicity 
due to interchange instabilities operating in the regime 
where ~ rc cannot be ruled out as a mechanism for 
the 1 Hz QPO. 

4.2.3. Magnetic reconnection instabilities 

As matter moves within the disk (radially and az- 
imuthally) it can drag magnetic field lines along with 
it. The sheared field lines can temporarily impede accre- 
tion until reconnection establishes a normal fiow again. 
The resulting quasi-periodic accretion flow would lead to 
a corresponding quasi-periodicity in the lightcurve, pro- 
vided that the accretion funnel and hot spot can respond 
on a 1 s time scale. 

Magnetic reconnection is one of the most plausible 
mechanisms for Type II bursts (cf. § 4). Type II bursts 
are thought to occur in systems where magnetic inhi- 
bition causes accreting matter to build up in a reser- 
voir outside the magnetosphere. Once a sufficient over- 
density of material has accumulated, instabilities cause a 
catastrophic breach of the magnetospheric hammock, re- 
sulting in sudden bursts of accretion (Lewin et al. 1976). 
Dramatic changes in QPO properties immediately before 
Type II bursts, with no detectable change in the spec- 
trum, have been interpreted as indicating that the QPOs 
in the persistent emission are generated within the fuel 
reservoir (Dotani et al. 1990). To what extent this phe- 
nomenon is relevant to our 1 Hz QPO is unclear. In the 
following discussion we explore several mechanisms that 
can produce 1 Hz oscillations. 

Aly & Kuijpers (1990) used reconnection to explain 



the QPOs observed in the Rapid Burster, and discussed 
how the disk would be broken up in 'blobs' by magnetic 
reconnection instabilities. They predict a frequency a 
few times the beat frequency between the rotation rate 
of the star Vs and the Keplerian frequency at the inner 
edge of the disk uk- This implies that the frequency of 
the QPOs should pass through zero at the point where 

An argument in favor of this mechanism is that has 
been observed in simulations (Goodson, Winglee & 
Bohm 1997; Goodson, Bohm & Winglee 1999; Good- 
son & Winglee 1999, and Romanova et al. 2005). Ro- 
manova et al. (2005) observed a strong quasi-periodicity 
associated with outflows of matter when the system was 
in the propeller regime. Unfortunately Romanova et 
al. (2005) do not model the emission mechanisms, so 
whether outflows or discrete accretion woiild genuinely 
lead to high amplitude QPOs in the X-ray emission is 
not clear. Ustyugova et al. (2006), using MHD simula- 
tions, showed that rapidly rotating accreting stars have 
strong periodicity of this type linked to strong outflows. 
Observationally, there is evidence for jet formation from 
the radio detections in the decay of the 1998 outburst 
and the peak of the 2002 outburst (Gaensler et al. 1999, 
Rupen et al. 2002), although the latter cannot be related 
with the propeller onset since it occurs at the maximum 
fluxes observed for .11808. 

The frequency, for this mechanism, should depend on 
accretion rate (as observed. Fig. 7). It should however 
occur at all accretion rates: current models suggest no 
means of conflning this mechanism to a narrow range of 
accretion rates as we observe for the 1 Hz modulation 
(§ 3.3). Another argument against this mechanism is 
that this type of instability should occur independently 
of field misalignment angle, so should be seen in all the 
AMXPs and the non-pulsating LMXBs (assuming that 
they all have an external magnetic field). 

4.2.4. Thermal/viscous/radiation-driven instabilities 

The inncT rc\a;ion of an a disk (Shakura & Syunyaev 
1973) is unstable to thermal and surface density pertur- 
bations (Pringle, Rees & Pacholczyk 1973; Lightman & 
Eardley 1974). Both types of instability can arise when 
radiation provides the major contribution to the total 
pressure (Shakura & Sunyaev 1976). Frank, King & 
Raine (2002) derive the various time scales in operation: 
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The time scales are defined as follows: is the dynam- 
ical time scale in the disk; Tz is the time scale on which 
deviations from hydrostatic equilibrium in the z-direction 
get smoothed out; Tth is the thermal time scale, that is 
the time scale for re-adjustment to thermal equilibrium, 
if, say, the dissipation rate is altered; and Tyisc is the 
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time scale on which matter diffuses through the disk due 
to the effect of viscous torques. Note that the canonical 
value suggested by numerical simulations for a is 0.01. 
Thus the dynamical time scale is of the order 10~^ s, the 
thermal time scale is 0.01 s and the viscous time scale is 
of the order 100 s assuming that the appropriate length 
scale is comparable to the inner disk radius R ^ Tc ^ 10 
km. A study by King, Pringle & Livio (2007) shows 
that for many LMXBs a ~ 0.1, an order of magnitude 
larger than the value suggested by numerical simulations 
of disks. Therefore the observed time scales are expected 
to be shorter than those calculated from numerical simu- 
lations. If the instability were confined to a narrow inner 
annulus of the disk, then this might be consistent with 
a 1 Hz frequency, since the viscous time scale falls for 
shorter lengthscales. In addition, if the instability trig- 
gers at the onset of the propeller regime when the inner 
edge of the disk puffs up to sustain accretion, this might 
explain the rarity of the phenomenon. 

Instabilities which are known to operate at high accre- 
tion rates, such as thermal-viscous instabilities in radia- 
tion dominated disks (Taam & Lin 1984) and radiation- 
driven instabilities (Fortner, Lamb & Miller 1989; Miller 
& Park 1995), can be immediately ruled out since the 1 
Hz QPO is observed at mass accretion rates of 0.001 — 
0.01% Eddington. 

An instability that might be relevant is the ionization 
instability that is thought to put the system into out- 
burst in the first place (see Lasota 2001 for a detailed 
review). We know that the 1 Hz QPO appears at the 
end of the outburst, where current accretion disk models 
predict a transition from the hot to the cold state. The 
ionization instability might enter a marginally stable os- 
cillatory state when the disk is on the verge of flipping 
between hot and cold regimes at luminosities of ^1% 
Eddington. The question is then why it would not oc- 
cur always when the source is in the required luminosity 
range, and why not also in other transient LMXBs. Fine 
tuning by requiring this marginal state to coincide with 
for example the onset of the propeller regime and the as- 
sociated changes in the disk structure as discussed above 
might resolve this. Therefore, the ionization instability, 
although unlikely, cannot be ruled out. 

4.2.5. Spruit-Taam instability 

The Spruit-Taam instability involves radial perturba- 
tions at the magnetospheric boundary (the inner edge of 
the accretion disk according to Spruit & Taam 1993). It 
relies on the viscosity in the disk to work, and does not 
require any shearing of the magnetic field. As discussed 
in Section 4.1, once Vyti accretion is only possible if 

the density at the inner edge of the disk rises. A small 
perturbation of the the disk radius away from — 
will be immediately damped and the disk radius will re- 
turn to the "equilibrium" position = rc where inner 
disk edge and magnetosphere have the same angular ve- 
locity. 

However, in the early propeller regime there exists a 
marginal state. When the inner accretion disk is at equi- 
librium, a given corresponds to a specific density. If 
Tm moves inward a little bit, the boundary layer is 'over- 
dense' compared to the equilibrium for that smaller rm-, 
and so quickly empties out. The rapid flow of matter 
from the inner edge of the disk causes the rest of the 



co-rotating transition zone to empty out too. The rise 
in local m pushes in still further, reinforcing the per- 
turbation. Eventually however the innermost layers are 
devoid of matter, since the viscous time scale further 
out in the disk is too slow to have replenished the inner 
regions. At this stage the boundary has to move back 
outward. Matter then accumulates again until the cycle 
restarts. 

The time scale on which this instability operates is re- 
lated to the viscous time scale just outside the corotating 
transition zone of the disk, since this sets the time scale 
for replenishment of the reservoir. However, this time 
scale is not related in a simple way to the viscous time 
scale as it also depends on parameters like the average 
accretion rate and the steepness of the transition between 
disk and magnetosphere. 

According to eq.(6), the viscous time scale at rc for 
J1808 is ~ 100 s, 2 orders of magnitude too long to ex- 
plain the 1 Hz QPO. However this neglects the depen- 
dence of the time scale on the additional parameters and 
assumes that the appropriate length scale is compara- 
ble to the inner disk radius rc- If a shorter length scale 
were involved then a shorter viscous time scale would 
be possible (Spruit & Taam 1993). To obtain a 1 sec- 
ond time scale, the length scale of the region of activity 
would need to be ~ 1 km, comparable to the scale of the 
inner "puffed up" regions in the steady state disk models 
calculated by Rappaport, Fregeau & Spruit (2004). The 
advantage of this mechanism is that the oscillatory state 
is expected to be active only in a very narrow range of 
radii from ~ rc to = l-5rc and hence a small range 
of accretion rates, as we observed for the 1 Hz QPO (see 
§ 3.2). 

The Spruit-Taam instability could modulate the accre- 
tion flow at high amplitude, which would flt the obser- 
vations of very high fractional rms amplitudes for the 
QPO. The instability would also be compatible with the 
continued presence of accretion-powered pulsations, since 
accretion could still be funneled even if the inner edge of 
the disk were oscillating. Finally, the frequency of the 
instability has a weak dependence on the mass accretion 
rate (see Fig. 4 in Spruit & Taam 1993), rising or falling 
whether is greater or less than Tc- This weak depen- 
dence has been observed in J1808 (Fig. 7) with the fre- 
quency rising with X-ray flux, thus suggesting > rc. 

4.3. The mechanism for the 1 Hz QPO 

In summary, most of the mechanisms examined cannot, 
based on our current understanding of how they work, 
explain key features of the 1 Hz QPO (see Table 2). The 
mechanisms that remain plausible are all associated with, 
or fine-tuned by, the onset of the propeller regime. 

There are a number of other pieces of evidence (cf. § 1) 
that also point to major changes in the accretion envi- 
ronment at the luminosity where the 1 Hz QPO sets in 
(Wijnands et al. 2001, Wijnands 2003, Campana et al. 
2008) - changes which might be explained by the on- 
set of the propeller. In addition there are timing results 
suggesting a major change in disk structure aroiuicl this 
time, such as the ~ 0.2 phase drift in the fundamen- 
tal (arguing for a major change in the disk environment 
around this time), the change in the soft lag behavior 
(Hartman et al. 2009b), and the (debated) detection of 
an accretion torque (Burderi et al. 2006; Hartman et al. 
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TABLE 2 
Mechanisms for the 1 Hz QPO 



Model Preq. Amp. Amp. /Energy Flux threshold Preq./Flux. Propeller 



Surface oscillations P N Y N • • • N 

Marginally stable nuclear burning N'-- ■■• Y ■■• N 

Disk obscuration Y N P • • ■ • • • N 

Interchange instability (rm ^ Tc) • ■ ■ ■ ■ ■ ■ ■ ■ N ■ ■ ■ N 

Interchange instability (r™ ~ Tc) ■ ■ ■ ■ ■ ■ ■ ■ ■ Y • • • Y 

Magnetic reconnection instability P • ■ ■ ■ ■ ■ N Y N 

Thermal / viscous instability P-— ••■ N •■• P 

Ionization instability P—' ••■ Y ••■ P 

Radiation instability N—' ••■ N ■■■ N 

Spruit-Taam instabiUty P Y • • • Y Y Y 



Note. — The table compares observed properties of the 1 Hz QPO to the various models discussed in 
§4. Y/N (yes/no) indicates that the model can/cannot explain the property in J1808 (e.g. at the accretion 
rates inferred for this source). The symbol P (possible) indicates that the model might be able to explain the 
property if certain conditions arc met. An empty field indicates that the model makes no specific prediction for 
that property, and that further studies arc required. Columns: (1) model, (2) 1 Hz QPO frequency, (3) high 
fractional amplitude, (4) energy dependence of the QPO fractional amplitude, (5) appearance of the oscillation 
below a certain fiux threshold, (6) dependence of the centroid frequency on the X-ray fiux. The last column 
shows whether the onset of the propeller regime is necessary for the model to be able to explain the observed 
QPO properties. Only models without N in any column remain viable, pending further study. 



2008). 

The mechanism proposed by Spruit & Taam (1993) 
seems to be the most promising candidate to explain the 
1 Hz QPO, although the precise details of the time scales 
for this instability in the situation when funnel flows are 
relevant remain to be worked out. It has a precise onset 
point associated with the early propeller regime, should 
remain relatively stable in frequency as accretion rate 
varies slightly and is only expected in a narrow range of 
accretion rates. 

Other mechanisms may also play a role, perhaps in con- 
cert with the Spruit-Taam instability. In § 4.2.2 we men- 
tioned that new classes of interchange instabilities might 
operate near the propeller transition, perhaps leading to 
sporadic accretion. In § 4.2.4 we discussed the possibility 
of the ionization instability triggering on short length- 
scales in the inner regions of the disk once the source 
enters the propeller regime. This possibility is particu- 
larly plausible if the disk is already close to the transition 
from outburst to quiescence. The ionization instability 
might reinforce the Spruit-Taam instability mechanism, 
and could also fine-tune the onset conditions for the 1 Hz 
QPO (sec § 4.4). The munber of empty fields in Table 2 
reflects the scale of the modeling work required to resolve 
these questions. 

It is hard to understand why the 1 Hz QPO does not 
appear during the faint re-flares in the 2008 outburst. 
8 out of 57 observations were in the 2-15 mCrab range 
during the re-flaring state. The reason why the 1 Hz 
QPO is not observed in these 8 observations is an open 
problem. Although poorly constrained, the 1998 out- 
burst exhibited a similar behavior, and on several occa- 
sions during the 2000, 2002 and 2005 outbursts the 1 
Hz QPO also remained undetected even for fluxes in the 
2-15 mCrab range, with fractional rms amplitude upper 
limits of ^ 10%. Clearly the 1 Hz QPO mechanism is 
not always triggered even in the 2-15 mCrab range in 
J1808. 



4.4. J1808 and the other AMXPs 



In order to enter the propeller regime, J1808 needs to 
be at the point where r„j ^ re- Equating the crude 
expressions given in eq.(3) and (4), we obtain a relation 
between accretion rate, magnetic field and spin rate. Fig- 
ure 10 shows the conditions for propeller onset for par- 
ticular combination of these parameters. Clearly this is 
very approximate, since it is based on the simplest esti- 
mates of and Tc, and ignores dependencies on mass 
and radius, but sufiicient to understand whether the pro- 
peller scenario is a realistic possibility. 

We plot the mass accretion rate values of three well 
known AMXPs: XTE J1807-294 (spin frequency 190 
Hz), J1808 (401 Hz) and ICR J00291+294 (599 Hz). 
The first object was chosen because its spin frequency is 
one of the lowest known among AMXPs and its outburst 
spans a wide range of luminosities. IGR J00291-I-294 
was chosen because its neutron star has the highest spin 
frequency known among AMXPs. 

The mass accretion rates used in Fig. 10 are calcu- 
lated for X-ray fluxes in the 2-10 kcV energy band, by 
assuming a neutron star mass of 1.4 Mq and an efii- 
ciency of 10% for the conversion of rest mass energy of 
the accreted material into X-ray flux. Since those mass 
accretion rates do not refer to bolomctric fluxes, they 
have to be considered lower limits. We also marked the 
bolometric peak luminosity of each source (as reported 
in Gierlinski et al. 2002, Falanga et al. 2005a,b) for as- 
sumed distances of 8.5 kpc (IGR J00291+294-'' and XTE 
J1807-294) and 3.5 kpc (J1808). The very broad range of 
luminosities of J1808 are observed thanks to the deeper 
observations of Swift-XBH (§ 3.1, Campana et al. 2008) 
and XMM-Newton (Wijnands 2003). 

For all three sources, the conditions for propeller onset 
should be encountered if the field strength is ^ 10® G. 
For J1808, with a spin of 401 Hz and an accretion rate 
that runs from a few percent of Eddington at peak, down 
to less than 0.001 % in the dips between the re-flares, 
the system must always enter the propeller regime at 

3 we note that IGR J00291+294 is not located toward the Galac- 
tic Center, therefore the assumed distance is arbitrary, and is made 
for an easier compajrison with XTE J 1807-294 
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Flc;. 10. — The accretion rate at which the propeller sets in for particular combinations of magnetic field and spin rate for the three 
AMXPs (dashed curves). The vertical arrows show the range of accretion rates over which the AMXPs have been observed and refer to 
X-ray fluxes. They have to be considered lower limits on the true mass accretion rates which are set by bolomctric fluxes. The estimated 
bolornetric fluxes at the peak of the outburst (where the mass accretion should be the highest) are marked with a black square (sec text 
for details). The thick dashed vertical line in J1808 indicates the range of mass accretion rates at which the 1 Hz QPO has been observed. 



same accretion rate, while for B < 10 G the system 
will not enter the propeller regime in the observed range 
of mass accretion rates. The range of magnetic fields 
is {B ^ 0.4-1.5 X lO^G) as reported by Hartman et al. 
(2008, 2009a). The range of accretion rates for which the 
1 Hz QPO appears (inferred from the 2-15 mCrab X-ray 
flux, § 3.2) lies just below this range. This coincidence is 
quite impressive since the accretion rates are lower limits. 

By looking at Fig. 10 we can infer that the main reason 
why the 1 Hz QPO has been observed in J1808 and not 
in other AMXPs might be related with the proximity of 
J1808 (3.5 kpc) with respect to the other AMXPs (as- 
suming they are all located at a distance close to 8.5kpc). 
From a visual inspection of Fig. 10 the non-observation 
of the propeller in IGR J00291 could be problematic if 
the maximum accretion rate is ~ 10% Eddington and 
the magnetic field is ~ 10*G. However, many uncertain- 
ties can play a role here: uncertainities on the source 
distance, the conversion of X-ray luminosities into mass 
accretion rates, the precise condition for the onset of the 
propeller, not to mention all the uncertainties related 
with the definition of niagnetospheric radius. In this 
sense Fig. 10 has to be taken as a qualitative picture, use- 
ful to understand the underlying behavior of this sources, 
but with too many uncertainties remaining to draw ro- 
bust conclusions. 

However, if we suppose that the instability is triggered 
in IGR J00291+5934 at the same mass accretion rates as 
in J1808, then the expected luminosity would be below 
the detection threshold of RXTE. None of the known 
AMXPs has been extensively monitored at low flux levels 
by Swift-XRT or XMM-Newton, both of which would be 
able to (easily) test this scenario. 

If J 1808 enters the strong propeller regime during the 
re-flares, a large outflow of gas is expected. Hartman et 



al. (2008) and di Salvo et al. (2008) observed an anoma- 
lously large orbital period derivative. A possible explana- 
tion requires a mass loss from the system of ~ 1O~^M0 
yr-i (di Salvo et al. 2008 and Burderi et al. 2009 for 
a description of this scenario). The onset of a strong 
propeller with outflows of matter from the system can in 
principle play a role in this. If however, the real explana- 
tion is different (see Hartman et al. 2009a for a discussion 
of alternative possibilities to the strong outflow scenario) 
then the propeller might have only a minor role, if any, 
on the long term evolution of the orbital period. 

A strong propeller might in principle explain the long 
term spin down of the neutron star in J1808 as pro- 
posed by Hartman et al. (2008). By assuming a B 
fleld of ~ lO^G, and assuming ^ 2rc we can cal- 
culate the amount of mass that needs to be ejected 
in a strong propeller regime to produce the long term 
spin down observed by Hartman et al. (2009a) {v ~ 
-5.5 X 10-16 Hzs-i): 

Nsd = Nprop -^277 IP = y^GMrmMej (7) 

where Nsd is the spin-down torque, Nprop is the pro- 
peller torque (defined as in Bildsten 1998) and M^j is 
the amount of mass expelled during the strong propeller 
stage. The mass ejection required during quiescence is 
~ 8 X 10~^^ Mq yT~^ . This mass outflow is still insuf- 
ficient to explain the large orbital period derivative ob- 
served in J1808. Moreover, even if the disk material is 
expelled at r„j ^ 2rc, the viscous dissipation of grav- 
itational potential energy would yield a luminosity of 
« 10"^^ergs-i, which is too high when compared to the 
quiescence luminosity of J1808 (~ 5 x 10^^ ergs-^, Cam- 
pana et al. 2002, Heinke et al. 2009). 

We have shown that the 1 Hz modulation has an effect 
on the 401 Hz pulse formation. Magnetic channeling is 
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expected to be easier in the valleys than in the peaks of 
the 1 Hz modulation, where the accretion rate is lower. 
We observed an increase of the 401 Hz pulsed fraction in 
the valleys in agreement with this. We have also found 
possible links between the 1 Hz modulation amplitude 
and the 401 Hz pulse phase. If the 1 Hz QPO reflects a 
change in the accretion flow, the position of the magnetic 
funnel can change accordingly, thus affecting the 401 Hz 
pulse phase. A better understanding of the process that 
generates the 1 Hz modulation would be fundamental to 
clarify this magnetosphere/disk interaction. 

5. CONCLUSIONS 

We have performed the first complete study of the 1 Hz 

modulation, its relation to the 401 Hz pulsations and the 
re-flares of SAX J1808. 4-3658 as observed over 10 years. 
Several common features are observed in the 2000, 2002 
and 2005 outbursts while the 1998 and 2008 outbursts 
have different properties, the most remarkable one being 
the absence of a strong 1 Hz modulation. 

We focused on the origin of the 1 Hz oscillation that 
sometimes dominates the re-flare lightcurve and we found 
that all viable candidate mechanisms are connected with 
the onset of the propeller stage. The most promising 
model discussed is the Spruit-Taam instability which ex- 
plains the stable 1 Hz frequency, its high amplitude and 
the narrow flux range of its occurrence in a natural way. 



and is also compatible with the simultaneous presence of 
401 Hz pulsations. 

Many open issues remain. It is not clear yet why the 
pulse phase drifts by 0.2 cycles during the fast decay 
just before the onset of the 1 Hz modulation. It is also 
unclear whether the pulse phase jumps observed on time 
scales of hours or less are related with the amplitude of 
the 1 Hz. However, it is likely that this pulse behavior 
reflects changes in the accretion flow onto the surface 
triggered by flow changes associated with the 1 Hz QPO 
mechanism. 

The reason why 1 Hz modulation has not been ob- 
served in other AMXPs might be the larger distance of 
the other sources making them undetectable in the rel- 
evant luminosity range. Future monitoring of low level 
flux states of other AMXPs will be very important to 
extend our comprehension of the 1 Hz modulation and 
further probe the onset of the propeller stage. 
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